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Abstract 

The physics of a Tonks-Giradeau Gas in the presence of a local potential is studied. In or- 
der to evaluate the single particle density matrix (SPDM) of the many-body ground state, the 
Wiger-Jordan transformation is used. The eigenvector with the largest eigenvalue of the SPDM 
corresponds to the " Bose-Einstein Condensate" (BEC) State. We find that the "BEC" state density 
at the positon of the local potential decreases, as expected, in the case of a repulsive potential. For 
an attractive potential, it decreases or increases depending on the strength of the potential. The 
superfluidity of this system is investigated both numerically and perturbatively. An experimental 
method for detecting the effect of an impurity in a Tonks-Giradueau gas is discussed. 
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I. INTRODUCTION 



The effect of impurity and disorder on many body system has been one of the main 
themes of condensed matter physics[l]. It plays an important role in our understanding of 
phenomena such as superconductivity, superfluidity, and Kondo physics. The rapid progress 
in ultra-cold atom experiments provides a unique platform for the study of such many-body 
systems. In these experiments, the external potential, the dimensionality, and the interaction 
can be tuned using external fields. Fascinating many body phenomena have already been 
realized in experiments [2], including recent realization of a To nks- Girardeau Gas [3]. 

A natural extension of this line of research is the study of the effect of impurities in 
these many body systems, both theoretically [7], and experimentally [8]. Here we study 
the Tonks-Girardeau gas in the presence of a local potential. A Tonks-Giradeau gas is a 
quantum gas consisting of hard core bosons. This system can be mapped to a free fermion 
system, allowing for an exact solution. The study of the effect of an impurity on such an 
exactly solvable system serves two purposes: first, it shows, in its simplest form, the interplay 
between interaction, superfiuidity, and impurity. Second, the theoretical predictions can be 
compared readily with experimental observables. 



II. MODEL OF THE TONKS-GIRARDEAU GAS WITH AN IMPURITY 



The Hamiltonian we use to describe such a system is of the tight binding type 

H (^alui+i + al_^^a^^ + -uajoo (1) 

i 

where a] and are the bosonic operators that create and annihilate a particle at site i, 
respectively. They obey the commutation relation aj,aj — ^ij- The particle interaction 

is included by imposing of — (^af^ = 0. The quantity t, which is generally negative, is 
the hopping amplitude between nearest neighbors and u characterizes the local potential 
strength at site zero. The behavior of the many body system is determined by the ratio of 
u/t. When \u/t\ < 1, the kinetic parts dominate and the particles distribute uniformly to 
lower the kinetic energy. On the other hand, when \u/t\ > 1, the local potential becomes 
important. In this regime, local charge fiuctuations are suppressed and the coherence among 
particles is degraded. This degradation is a main interest of this paper. We numerically 



evaluate the SPDM of the many-body ground state, and compute the effect of a local 
potential (impurity) on the spectrum of the sSPDM, on the "BEC" wavefunction, and on 
the superfluidity of the many-body ground state. 

In order to solve for the ground state of the system, the Hamiltonian(l) can be mapped 
to an non-interacting fermion Hamiltonian via the Wigner- Jordan Transformation 

= n e'^'^-'^-Q, 4 = 4 n e-^'^"«"" (2) 

a=l a=l 

Here the cj and q are fermionic operators, which satisfy the anticommunication rela- 
tion |cj,c]| = 6ij. Under this transformation, it is straightforward to show that the 
Hamiltonian[l] can be written in the following form 

if = i ^ (cJq+1 + cJ_^iQ^ + uclco (3) 

i 

The ground state \Gf) of this Hamiltonian corresponds to the well-known Fermi sea, in 
which fermions fill the single particle levels up to the Fermi surface. The corresponding 
ground state \Gb) of the boson Hamiltonian(l) can be obtained from \Gf) by symmetrizing 
the corresponding many body fermionic wave function. 



III. CONDENSATE FRACTION AND CONDENSATE WAVEFUNCTION 



We seek the one particle density matrix [9], which is the expectation value of ajaj in the 
ground state. 

Pij = {Gb\ ajoj \Gb) 

= {Gf\ 4 n e-^^""^« if e^'^^«'="Cj- \Gf) (4) 

a=l a=l 

The diagonal part of the density matrix, pa, gives the density (rii) = {GbIclIch \Gb) 
at site i, while the off diagonal part gives the coherence in the many body ground state 
for different sites. In the uniform case, the coherence pij depends only on the difference, 
|i — j|. In the presence of a local potential, however, the coherence depends on both i and 
j separately. In addition to that, pij is much smaller than its uniform counterpart when 
i and j are on different sides of the local potential. Diagonalizing the density matrix gives 
a set of eigenvectors corresponding to the reduced one-particle states. The corresponding 
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eigenvalues represent the occupation probability for the corresponding reduced one particle 
states. The one with a significant large eigenvalue is the "BEC" State [9] (Note that in 
ID, there should be no Bose Einstein Condensation in the thermodynamical limit. In the 
particular case of the Tonks-Girardeau case, the particles that condense into a single state 
are calculated to be order of \fN [4]. However, this state, compared with other states, is 
the only state that is significantly occupied. With this clarification, we will call that state 
the BEC state and the corresponding occupation the condensate fraction) . 

The evaluation of this density matrix is not simple even when the ground state is known. 
Here we adopt the technique developed by M. Rigol and A. Muramatsu [6]. This method 
allows us to calculate a system with up to a thousand lattice sites. We focus on the low 
occupation limit only, which corresponds to a continuum case. We have verified that 9 
particles in 100 lattice sites already converges very well to the continuum case. In the 
following calculation, unless otherwise stated, all the calculations are based on a system of 
9 particles in one dimensional lattice of 100 sites, with an impurity located at site zero. 
Periodic boundary condition is imposed on this one dimensional chain. 

Since all the physics is determined by the ratio u/t, in the following numerical calculations, 
we fix t = —1 and vary u. Fig [la] shows single particle density matrix spectrum for difi^erent 
u's.The significant peak corresponds to the BEC state occupation or condensate fraction. 
The introduction of an impurity lowers the condensate fraction, as shown in Fig [lb]. The 
attractive and repulsive potentials almost equally deplete the condensate fraction for small 
impurity potential with the attractive potential having a slightly larger effect. The physics 
that gives rise to this behavior, as is shown in the following, is that an attractive local 
potential has a stronger effect on decreasing the coherence between particles. 

Since only the BEC state is significantly occupied, the BEC state determines the most 
important features of the many-body system. It would be useful to look at the BEC wave- 
function itself (see figure [2a]). For a repulsive local potential, we find that the BEC density 
decreases near the impurity. In the case of an attractive potential, for m > — 1, there is an 
increase in the probability of the BEC density at the impurity site. For the case of -u < — 1, 
in contrast to what one might expect, there is a decrease of the BEC density at the impu- 
rity site. This feature actually arises from the competition between two effects: the single 
particle effect, i.e. the potential attracts particles, and the many body effect, i.e. the im- 
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FIG. 1: (a) shows the SPDM spectrum at different local potential strength u. Here we have 
normalized the population to be one. One of the reduced one-particle states is overwhelming 
occupied compareing with all the others, and it is identified as the BEC state. We show only part 
of spectrum and omitting parts which are negligibly small, (b) shows the condensate fraction as a 
function of the local potential strength. The attractive potential is seen to have a slightly larger 
effect in decrease the condensate fraction. 

purity decreases the coherence among particles. We have also included the particle density 
plot in Fig [2b]. It shows an increase of particle density for any attractive potential. For 
the BEC density, we see that, for each attractive potential, there is a peak at the impurity 
site corresponding to the local bound state. However, due to the lack of coherence between 
bound and extended states, the bound state is less likely to participate the BEC state, which 
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results in a overall decrease of the BEC density at site zero compared with the uniform case. 
This is the main observation of this paper. To be more specific on how the impurity de- 
creases the coherence among particles, we plot some of the relevant off-diagonal elements of 
the single particle density matrix (see Fig [3]) In particular, we take the coherence between 
an arbitrary site and the fifth site next to the local potential site zero as an example. We 
find that, given the same distance, the coherence between particles on different sides of the 
impurity is much smaller than that of the particles located on the same side of the impurity. 
For the same magnitude of the attractive and repulsive potential, we find that the effect of 
the attractive and repulsive potential are roughly the same, with the attractive one having 
a stronger effect in decreasing the overall coherence among particles. This actually explains 
the fact the attractive potential has a stronger effect on decreasing the condensate fraction. 
The only place that the attractive potential resulting in a large coherence than the repulsive 
one with the same magnitude is near the impurity site. This is due to the presence of a local 
bound state, which effectively increase the probability of finding particles at the impurity 
site. This bound state is well known [10] in the tight binding Hamiltonian. Its eigenenergy 
is = —\/v? + and its wavefunction localizes in space as exp [—a (E) \n\] with n being 
the site number and 



corresponds to the inverse of the characteristic length. The nice overlap of the coherence 
peak and the bound state wavefunction verifies our argument for the increase of coherence 
in the vicinity of impurity. The small peak in the impurity site is purely a single particle 
effect. 

With the picture of the impurity introducing decoherence, one would suggest that in 
higher dimensions, since the impurity has a weaker effect in decreasing coherence among 
particles, a smaller effect of the impurity should be found. We have verified this point 
numerically by extending our calculation to two dimensions The condensate fraction is shown 
in Fig [4]. At dimensions greater than one, there is no simple mapping from bosons to 
fermions. An exact diagonalization has been done for this case which limits the calculation 
to a system of five particles in a three by three lattice. 
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FIG. 2: (a) The BEG wavefunctions for different strengthes of the local potential are shown. In 
the special case of u=0, the BEG wavefucntion is constant all over the lattice. Note that the BEG 
wavfunction corresponding to u = —1 has a lower value at the impurity site comparing with the 
u = uniform case, (b) The particle density is shown with respect to lattice sites for different u's. 
Here we see for attractive u, the particle density at the impurity site is always larger than that of 
the uniform case. 

IV. SUPERFLUIDITY 

It is well known that BEC is neither necessary nor sufficient for the existence of a super- 
fluid. As mentioned by Lieb and Sieringer[5], a particular example is the Tonks- Girardeau 
gas. Even it does not Bose condense, the system without impurity actually exhibits super- 
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Coherence between different sites 




FIG. 3: The off diagnal SPDM elements tliat measure the coherence between lattice sites -5 and any 
arbitary lattice sites is plotted for u=-l and u=l, respectively. The coherence with the negative 
sites are relatively larger than the coherence with the positive sites. The smaller peak of the 
coherence for the attractive potential between any sites and the sites zero can be explained by the 
increase of the density at impurity site, due to the presence of a bound state. The bound state 
wavefunction is plotted and it conincides well with the peak of the coherence in the impurity site. 



fluidity. It is interesting to investigate how the superfluidity is decreased in the presence of 
an impurity. We emphasize that we use the word superfluidity strictly in the sense of the 
following phenomenological definition 

E{v) E{0) 1 



N N 2 



fsmv'^ + O (v*) 
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Condensate Fraction in Two Dimension 
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FIG. 4: This figure shows the condensate fraction in a three by three lattice with an impurity 
located in the center. The impurity is charaterized by a local potential strength u. The total single 
particle density spectrum is normalized to one. We see that in two dimension, the condensate frac- 
tion is larger than its one dimensional counterpart. Morever, the relative change of the condensate 
fraction due to the presence of the local potential is much smaller than the one dimensional case. 

where E (v) is the ground state energy of the many body system in the presence of a 
perturbative velocity field v, fg is the superfluid fraction, and m is the mass of the particle. 
The velocity field is introduced by imposing a twisted boundary condition [5], which amounts 
to a phase jump e*'^, = whenever the wave function passes through the boundary. We 
restrict (f < n to yield a single valued function. Note that since the definition is based only 
on the static properties of the many body system, it can tell us only whether the ground 
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state has the property of superfluidity, and it cannot predict the stabihty of the superfluidity. 

The result of numerical calculations is shown in Fig [5] . We see that without impurity, 
the system exhibits superfluidity with superfluid fraction close to unity. The degree of 
degradation on the superfluid fraction produced by impurity depends on both the local 
impurity strength and the size of the system. In a large system, the superfluidity drops 
sharply with the presence of the local impurity. Note that as far as the spectrum is 
concerned, the fermi ground state is equivalent to the boson ground state. This allows an 
understanding of the superfluidity in a single particle picture. As the lattice size increase, 
the particle kinetic energy goes hke j2 while the local potential energy goes hke j^. Therefore, 
in the limiting case of large L, the local potential always dominates the kinetic energy. For a 
repulsive potential, the particle has to hop through the barrier. For an attractive potential, 
since the local bound state is always occupied, for hard core bosons, in order to get through 
this potential, a potential barrier of order is also present, which gives rise to the sharp 
drop of superfluidity in the presence of a local potential. 

This problem can be understood quantitatively by considering the following single particle 
model H = + 9^ {x)-, with x defined in the regime [— L, L] with the twisted bound- 

ary condition. As an example, we consider only two states with energy close to ^^^2 ■ To 
calculate the eigenenergy with different boundary condition is elementary. It is just the one 
dimensional piecewise potential problem. The modification of eigenenergy due to the pres- 
ence of the velocity field is [E (v) - E (0)] /2 = \mv'^ (1 - O.OOSc/^L^) [11] in the perturbative 
hmit that gL is small. This result shows that the superfluid fraction is /s = 1 — O.OOSg'^L^. 
There is no modification to first order in gL, which is consistent with discussions concerning 
repulsive and attractive potentials. 

V. DISCUSSION AND CONCLUSION 

The effect of an impurity on the Tonks-Girardeau Gas can be studied in an experiment 
using an optical lattice. The local potential we model here may be realized in an experiment 
by introducing a laser beam, an impurity atom, or an ion at a particular lattice site. The 

real space density distribution will be the same as the free fermion in the presence of a 
local potential, which contains only information on the diagonal part of the density matrix. 
The non-trivial part of the density matrix, the off diagonal coherence, can be detected by 



10 



Superfluidity witli different u's 
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FIG. 5: Superfluid fraction is plotted as a function of different local potential strength u. The case 
of nine particle occupying 100 and 200 lattice sites are shown. The larger the system, the sharper 
the superfuid peak is found 

measuring the momentum distribution, defined as 

N/2 
m,n=-N/2 

We have ignored the locahzed Wannier function profile. 

In figure [6] we show peak value density in momentum space for different impurity 
strengths. The highest momentum density peak corresponds to the case with no impurity. 
The attractive potential has a larger effect on the broadening compared with the repulsive 
one. The superfluidity of the system can be measured by imposing a velocity field on the 
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Momentum Peak Values For Different u's 
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FIG. 6: The peak value of the density distribution in the momentum space is shown. The maximum 
of zero momentum occupation is reached in the uniform case. 



many-body system. This can be realized in the experiment by changing of the external 
potential with time. The measurement can be done by looking at the damping motion of 
the particles. 

In summary, we have presented a numerical calculation of a Tonks-Girardeau Gas in the 
presence of a local potential. It is shown that the impurity decreases the occupation of 
condensate fraction. The BEG density at the local potential decreases in a strong attractive 
potential, as oppose to particle density, which always increases for attractive potentials. 
The superfluidity is also degraded by the impurity, the degradation scales like g'^L'^ in the 
small gL limit. This effect of impurity on a many-body system can be measured by looking 
at the momentum distribution and the system response an external velocity field. 
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